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PATENT 

Attorney Docket No.: 040092-0275 10US 
Client Ref. No.: MS-01437 and MS-01439 

METHOD AND SYSTEM FOR MEASURING DIFFERENTIAL 
SCATTERING OF LIGHT OFF OF SAMPLE SURFACES 

CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001 ] This application claims priority to U.S . Provisional Patent Application 60/429, 181, 
filed November 25, 2002, titled METHOD AND SYSTEM FOR MEASURING 
DIFFERENTIAL SCATTERING OF LIGHT OFF OF SAMPLE SURFACES, and is 
incorporated herein for all purposes. 

BACKGROUND OF THE INVENTION 
[0002] The present invention relates generally to an optical system, and more particularly 
related to a method and system for predicting the directions and intensities at which light 
scatters from a surface, such that a scattering prediction includes a substantially continuous 
solution over a set of scattering angles that includes specular and non-specular scattering 
angles, wherein the scattering prediction may be used to generate improved optical systems 
and improved scattering simulations. 

[0003] Scattering of stray light in optical devices, such as lenses, telescopes and the like, is 
a phenomenon that optical designers generally prefer to reduce. Stray light scattering 
generally refers to the scatter of stray light in an optical system from a surface, such as body 
portions of lens or telescope. Stray light in an optical system includes light that has separated 
from an optical signal and often travels along undesired paths. For example, stray light in the 
lens of a camera can cause photographic images to be formed that have aberrant light features 
that were not present in a photographed scene. 

[0004] Designers of optical system typically aim to reduce stray light generation and, once 
generated, to lower its degrading effect. To lower the degrading effects of stray light scatter, 
optical designers strive to generate relatively accurate predictive models of light scattering 
from surfaces. In generating relatively accurate models of light scattering from surfaces, 
stray light scattering may be better understood and compensated for in the design of optical 
systems. Models for stray light scattering can be incorporated in optical design programs that 
show the deleterious effects of stray light scattering as predicted by the models. Stray light 
scattering estimations calculated prior to building an optical system, can be used to improve 
final designs to build improved optical systems. 



[0005] While understanding stray light scattering in optical systems may be used to 
improve optical systems through improved optical system designs, models for stray light 
scattering can also be used in light simulation systems, such as in a computer graphics 
program (e.g., video games, and movie animations). Relatively accurate prediction of stray 
5 light scattering from a surface can be used in a computer graphic to simulate scattering from a 
surface that is substantially realistic. 

[0006] Predictive light scattering has historically been a diagnostic tool in the field of 
surface roughness measurements. It was realized relatively early in the study of scattering 
that the distribution of scattering from a surface can be related to the power spectral density 
10 of variations in surface heights of the surface. Generating relatively accurate models of 

power spectral density as related to surface height variation has provided for the extraction of 
surface height information from measured power spectral densities to characterize surfaces. 
See, for example, the studies of James Elmer Harvey reported in Light-Scattering 
Characteristics of Optical Surfaces, dissertation submitted to the Faculty of the Committee on 
15 Optical Sciences, The University of Arizona, 1976, herein Harvey '76. While predictive 

models of scattering have been used with some success in predicting scattering from surfaces 
and thus providing extraction of surface height information (e.g., degree that a surface 
resembles a diffraction grating pattern) from a surface, predictive models of scattering have 
generally provided non-continuous solutions, for example, between specular and non- 
20 specular regions of scatter. Scattering regions for which known models of scattering break 
down are often compensated for with empirical scattering data or intuitively estimated. 
However, compensating empirically or intuitively for the shortcomings of known models of 
scattering tends to be relatively time consuming and often costly as technicians and/or 
engineers often fill in manually, regions for which known scattering models break down. 
25 [0007] Measurement techniques and models for stray light scattering have historically been 
non-differential techniques. Non-differential techniques and models yield results that often 
include scattering information of a particular lab setup as well as scattering information of a 
sample. Separation of lab setup information from sample information in both measurement 
techniques and models of stray light scattering is often intuitively or empirically carried out. 
30 Addition steps of intuitive and/or empirical processing for isolating scattering information 
from lab setup information is often costly, time consuming, and may produce less than 
optimal results. 

[0008] Accordingly, what is needed are improved methods and systems for predicting 
scattering of stray light from surfaces and more specifically what is needed are models of 
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stray light scattering that generate continuous scattering solutions over a set of scattering 
angles that includes specular and non-specular scattering angles, wherein the scattering 
predictions may be used in computer design and computer simulations to generate improved 
optical systems and improved scattering simulations. 

BRIEF SUMMARY OF THE INVENTION 
[0009] The present invention provides an optical system, and more particularly provides a 
method and system for predicting the directions and intensities at which light scatters from a 
surface, such that a scattering prediction includes a substantially continuous solution over a 
set of scattering angles that includes specular and non-specular scattering angles, wherein the 
scattering prediction may be used to generate improved optical systems and improved 
scattering simulations. 

10010] According to one embodiment, a computerized method is provided for estimating 
scattering of electromagnetic radiation from a surface. The method includes providing a 
distribution expression that includes a first integral over a source solid angle, a second 
integral over a sample area, a third integral over detector solid angle, and an integrand that 
includes a differential-scattering profile; approximating the first and second integrals to be 
the second integral, wherein the source electromagnetic radiation is approximated to be 
collimated; approximating the second and third integral to be the third integral, wherein a 
detector for detecting the electromagnetic radiation scattered from the surface is 
approximated to be a point detector; transforming the coordinates of the third integral over 
detector solid angle to first and second dimensions in cosine space to form a fourth integral, 
wherein the surface is approximated to be shift invariant; integrating over the first dimension 
of the fourth integral; differentiating the fourth integral with respect to the second dimension 
to generate the differential-scattering profile; and generating an optical system design based 
on the differential-scattering profile. According to a specific embodiment, the distribution 
expression includes a bidirectional reflectance distribution function (BRDF). According to 
another specific embodiment, the differential-scattering profile is a continuous solution 
representing an algebraic model of specular scattering and non-specular scattering of the 
electromagnetic radiation from the surface. According to another specific embodiment, the 
differentiating step includes deconvolving the fourth integral. According to another specific 
embodiment, the step of approximating the first and second integrals to be the second integral 
includes approximating a one-to-one correspondence between a differential element of the 
source electromagnetic radiation and a differential surface area of the surface. According to 



another specific embodiment, the step of approximating the second and third integral to be 
the third integral includes approximating that electromagnetic scattered from a differential 
surface area sources is incident on the point detector. 

[0011] According to another embodiment, a computerized method is provided for 
estimating scattering of electromagnetic radiation from a surface. The method includes 
providing a distribution expression that includes a first integral over a source solid angle, a 
second integral over a sample area, a third integral over detector solid angle, and an integrand 
that includes a differential-scattering profile; approximating the first and second integrals to 
be the second integral, wherein source electromagnetic radiation is approximated to be 
collimated; approximating third integral to be one based on detecting the electromagnetic 
radiation scattered from the surface at an imaging detector; transforming the coordinates of 
the second integral over the sample area to first and second dimensions in cosine space to 
form a fourth integral, wherein the surface is approximated to be shift invariant; integrating 
overthe first dimension of the fourth integral; differentiating the fourth integral with respect 
to the second dimension to generate the differential-scattering profile; and generating an 
optical system design based on the differential-scattering profile. 

[0012] According to another embodiment, an optical system is provided and includes a 
collimated beam of electromagnetic radiation configured to illuminate a sample surface, the 
sample surface being shift invariant; an imaging detector configured to collect 
electromagnetic radiation scattered from the sample surface, the imaging detector configured 
to collect the scattered electromagnetic radiation at a plurality of scattering angels to generate 
a scattering profile; and a computer device configured to generate an estimated-differential- 
scattering profile and compare the scattering profile and the estimated-differential-scattering 
profile to generate an optical system design, wherein the estimated-differential-scattering 
profile is a continuous solution of an differential model of spectral scattering and non-spectral 
scattering derived from a deconvolution of a bidirectional reflectance distribution function 
(BRDF). According to a specific embodiment, an expression for the BRDF includes a first 
integral over a source solid angle, a second integral over the sample surface, a third integral 
over detector solid angle, and an integrand that includes the estimated-differential-scattering 
profile; the first and second integrals are approximated to be the second integral based on the 
source electromagnetic radiation being in the form of the collimated beam; the third integral 
is approximated to be one based on the detector being an imaging detector; the second 
integral is are transformed from an integral over detector solid angle to a fourth integral over 
first and second dimensions in cosine space based on the sample surface being shift invariant; 



and the fourth integral is integrated with respect to the first dimension and deconvolved with 
respect to the second dimension to generate the estimated differential-scattering profile. 
[0013] Numerous benefits may be achieved using the present invention over conventional 
techniques. For example, the invention provides for estimating differential-scattering profiles 
5 from surfaces, wherein the differential-scattering profiles are continuous solutions between 
specular and non-specular regions of scattered electromagnetic radiation scattered from a 
surface. Continuous differential-scattering profiles may be used to build improved optical 
systems, such as lenses and telescopes, that have reduced scattering or for which scattered 
radiation can be compensated for to lower the deleterious affects of the scattered radiation. 

1 0 Continuous differential-scattering profiles may also be used to build improved computer 
graphics programs (e.g., computer game graphics) to more accurately represent 
electromagnetic radiation reflected form graphical computer objects. Depending on the 
specific embodiment, there can be one or more of these benefits. These and other benefits 
can be found throughout the present specification and more particularly below. 

15 [0014] A further understanding of embodiments of the invention will be understood by 
reference to the following detailed descriptions and claims, and to the appended drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0015] FIG. 1 is a block diagram of a computer system in which embodiments of the 
20 present invention may be implemented; 

[0016] FIG. 2 shows a sample having a surface from which incident electromagnetic 
radiation may be scattered; 

[0017] FIG. 3 is a high-level flow chart having steps for deconvolving the bidirectional 
reflectance distribution function (BRDF) according to an embodiment of the present 
25 invention; 

[0018] FIG. 4 shows electromagnetic radiation emanating form a unique differential solid 
angle and traveling essentially parallel onto a differential element of the surface; 
[0019] FIG. 5 shows a scattering distribution in direction-cosine space for a shift-invariant, 
isotropic, and homogenous surface; 
30 [0020] FIG. 6 shows the reflection coordinate system (r, (f> ) for scattered electromagnetic 
radiation scattered from an illumination spot that is illuminated by incoming electromagnetic 
radiation; 
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[0021] FIG. 7 is a high-level flow chart having steps for deconvolving a generalized BRDF 
expression according to another embodiment of the present invention; 
[0022] FIG. 8 is a high-level flow chart having steps for deconvolving BRDF 
measurements to generate a distribution-scattering profile (DSP) according to an embodiment 
5 of the present invention; 

[0023] FIG. 9 is a high-level flow chart having steps for deconvolving BRDF 
measurements to generate a DSP according to another embodiment of the present invention; 
[0024] FIG. 10 shows an optical system in which embodiments of the invention may be 
implemented; 

1 0 [0025] FIG. 1 1 shows an example of a DSP that may be calculated from interpolated BRDF 
data according to an embodiment of the present invention; and 
[0026] FIG. 12 shows another example of a DSP that may be calculated from non- 
interpolated BRDF data according to another embodiment of the present invention. 

1 5 DESCRIPTION OF SPECIFIC EMBODIMENTS 

[0027] The present invention provides an optical system, and more particularly provides a 
method and system for predicting the directions and intensities at which light scatters from a 
surface, such that a scattering prediction includes a substantially continuous solution over a 
set of scattering angles that includes specular and non-specular scattering angles, wherein the 

20 scattering prediction may be used to generate improved optical systems and improved 
scattering simulations. 

[0028] FIG. 1 is a block diagram of a computer system 100 in which embodiments of the 
present invention may be implemented. A specific embodiment of the invention is 
implemented on a computer system 100 having a processor 105, a system memory 110, 

25 output devices 115, input devices 120, a disk memory 125, I/O (input output) ports 130 and 
an interconnecting device 135. Processor 105 may be implemented in a variety of formats, 
such as, but not limited to, a microprocessor, a microcontroller, a microcomputer, embedded 
logic or other processor types. Processor 105 may be a microprocessor manufactured, for 
example, by Intel Corporation, Motorola, Inc., or Advanced Micro Devices, Inc. System 

30 memory 1 10 may include EPROMs, EEPROMS, flash memory, SRAMs, DRAMs, cache 
memory and the like. Output devices 115 may include a variety of device types, such as 
CRTs, liquid-crystal display panels, printers, computer networks, an audio playback device 
and the like. Input devices 120 may include a variety of device types, such as a microphone, 
a keyboard, a computer network and the like. I/O ports 130 may include serial and parallel 
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ports of various designs, such as skuzzy ports, universal serial buses (USB) and the like. It 
should be understood that computer 100 may be local computer system, such as a personal 
computer, a distributed system with various modules located at a variety of locations 
connected by interconnect device 135 that may be a network, such as a WAN, a LAN, the 
5 Internet, or a PC system bus. A computer software program stored on system memory 110 
and/or disk memory 125 is configured to generate machine code instructions that implement 
embodiments of the present invention. While various aspects of the invention may be 
implemented in software, aspects of the invention may be implemented in hardware, 
firmware, control logic, or a combination of the foregoing. 

10 [0029] FIG. 2 shows a sample 200 having a surface 202 from which incident 

electromagnetic radiation 205 is scattered. Scattered electromagnetic radiation is identified 
by reference numeral 215. Incident electromagnetic radiation 215 is shown to be incident on 
surface 202 at incident angels cp* and Q\ (identified by reference numerals 220a and 220b, 
respectively) and scattered into scattering angles cpa and 0d (identified by reference numerals 

15 225a and 225b, respectively). Distribution of electromagnetic radiation into angles 9d and 9 d 
may be described by the bidirectional reflectance distribution function (BRDF). The BRDF 
is a generalized formula for modeling scattering, such as stray light scattering, of light from a 
source scattered from a surface. The BRDF may be written as follows: 

BRDF = - — — f f f AifiMjo^. 

p. Q hi JArea hd dQdA dCl d ' ' 

20 In the BRDF, Qi is the solid angle of incident electromagnetic radiation 205 incident on 

surface 202. Pi is incident electromagnetic radiation 205 and may represent incident power, 
incident intensity or the like. Qd is the solid angle of scattered electromagnetic radiation 215 
scattered from surface 202. Pd is the scattered electromagnetic radiation 215 and may 
represent scattered power, scattered intensity or the like. A is the area on which incident 

</p d (Q.,Q.,A) 

25 electromagnetic radiation 205 strikes surface 202. — — * s the differential- 

scattering profile (DSP) of scattered electromagnetic radiation 215. The BRDF is normalized 

by the factors — and — !— . 

Pi Cli 

[0030] Each of the integrals shown in the BRDF is a combination of two integrals. The 
integral over the differential solid angle d£l\ represents integrals over the incident angles q>\ 
30 and 0i, having differentials d<p\ and d&\ 9 respectively. The integral over the differential solid 
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angle dCld represents integrals over the incident angles (pa and 0a, having differentials dcp^ and 
J0d, respectively, and the integral over the differential area dA represents integrals over 
surface coordinates, such as x and y having differentials dx and dy, respectively. While dA is 
described as being a differential of x and y, which may be linear coordinates, dA may be a 
differential of polar coordinates or the like. 

[0031] According to one embodiment, the BRDF is simplified and deconvolved such that 

the differential-scattering profile ^ p d( Q " Q " A ) Q f sca ttered electromagnetic radiation 215 

dQ d 

is solved for. Deconvolution may be roughly thought of as inverting the integrals to extract 

dp d (n i9 n d9 A) 

the integrand (i.e., — ^— l - — d - — Simplifications of the BRDF include differential 

dil d 

simplifications, which yield a differential-scattering profile that is a continuous 
solution/model for scattered electromagnetic radiation 215 from specular regions into non- 
specular regions. Two exemplary methods for simplifying and deconvolving the BRDF are 
explained in detail below. Each exemplary method will be described in combination with 
respective high-level flow charts outlining steps included in the exemplary methods. 
[0032] FIG. 3 is a high-level flow chart having steps for simplifying and deconvolving the 
BRDF according to an embodiment of the present invention. It should be realized that the 
steps shown in FIG. 3 are not limiting on the invention as recited in the claims. Other 
techniques having fewer, substitute, and/or additional steps are within the purview of the 
invention and will be readily apparent to those of skill in the art. At 300, incident 
electromagnetic radiation is approximated as being essentially collimated. From a physical 
optics point of view, this means that the rays of incident electromagnetic radiation 205 are 
essentially parallel. In view of this approximation, rays of incident electromagnetic radiation 
emanating from a unique differential solid angle dCli are incident on a corresponding unique 
differential surface area dA. 

[0033] FIG. 4 shows rays 405 emanating from a unique differential solid angle 410 and 
traveling essentially parallel and incident on a differential element 420 of surface 202. As 
there is an essentially one-to-one relationship between unique differential source solid angles 
dQi and unique differential surface area dA, the contribution to the BRDF from these 
differentials is a delta function for each unique differential source solid angles dQi and 
corresponding unique differential surface area dA. Therefore, the integral over either the 
source solid angle or the integral over the surface area of the sample can be replaced with 
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one. In the embodiment presently described, the integral over the source solid angle will be 
set to one. An embodiment in which the integral over the surface area is set to one will be 
described below. Setting the integral over the differential solid angle dQ\ to one may be 
viewed as combining the integral of the differential solid angle </Qi and the differential 
5 surface area dA due to the interdependency of these differentials from the collimation 
assumption. 

[0034] As the integral over the differential solid angle </fii has been set to one, the 
expression for the BRDF has been reduced to a four dimensional integral that may be written 
as: 

10 BRDF = — — — f f ^^(^^>A) 

p. Q . JArea dA ^ * 

At 310, a point detector is assumed to collect scattered electromagnetic radiation 215 
scattered from surface 202. Detectors that may approximate point detectors include, for 
example, photomultiplier tubes, silicon diode detectors or the like. In view of this 
approximation, the integral over the differential scattering solid angle JQd becomes a delta 
1 5 function and may be set to one. As the integral of the differential scattering solid angle JQd 
is set to one, the BRDF may be rewritten as: 

BRDF = lJ-f ^ d M^^A) dA 

p. q. j a™ dA dn d 

At 320, sample surface 202 is assumed to be shift-invariant, isotropic, and homogenous. 
Shift-invariance implies that the shape of the differential-scattering profile does not changes 
20 in direction-cosine space with the angle of incidence and implies that 
</p d (Q /? Q rf ,A) 

^— - depends on \p -p Q \ . f3 is a coordinate in direction-cosine space and is a 

constant. Shift-invariance further implies that ^ p d( Q " Q <" A ) does not depend of surface 

d£l d 

position, surface orientation, or absolute scattering angle. \J3-J3 0 \ relates to scattering solid 
angle Qd via the relationship, |/?-y# 0 | = (sin 2 0i + sin 2 0d - 2sin9isindedcosA<|>d) ,/2 (this 
25 relationship is discussed in further detail below). 

[0035] In direction-cosine space, scattering intensities are equal on symmetric curves 
around the specular reflection. FIG. 5 shows a scattering distribution in direction-cosine 
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space for a shift-invariant, isotropic, and homogenous surface. Central circle 505 represents 
the specular reflection and each symmetric curve 510 around central circle 505 has an equal 
scattering intensity along its curve. At 330, the differential scattering distribution may be 
rewritten in direction cosines space as 

5 BRDF=— — — f fiMMi A . 

P. Q, jArea dA dQ. d 

[0036] In direction cosine space along radii that perpendicularly intersect central circle 505 
and symmetric curves 510, the intensity of the scattered electromagnetic radiation changes as 
a function of radius r. At 340, a coordinate transformation is made in direction cosine space. 
The coordinate system includes one coordinate <j) lying along the curves of constant \P~P 0 \ 
10 and the other coordinate r is perpendicular to symmetric curves 505 and 510 (i.e., in the 
|/?-/? 0 | direction). The variable r is a function of \J3-J3 0 \ . 

[0037] The integral along curves of constant |/?-/? 0 | are the path length of the curves (i.e., 
constants). A single integral is left over \/3-J3 Q \ that is simply an integral along the constant 
• curves of |/?-/? 0 1 . That is, dA may be written as d \ - fi 0 \ da , which yields the equation 

15 BRDF = — — f ^ dp * ^P~ Po ^ d\fi-P 0 \da. The integral along d\p-fi 0 \ arepath 
P. Q, jArea dA d£l d 

lengths of constant curves of |/?-y# 0 | . Both sides of the BRDF are integrated to yield an 

dp 1 dBDRF t 

expression for the differential-scattering profile: = - , wherein r is a 

r dQ I s ^(r)r dr 

coordinate that is a function of |/?-/? 0 | and £(r) is the range of the perpendicular coordinate. 
[0038] The following discussion describes a formalization for r and a formalization of a 

20 coordinate system. To build a coordinate system based on contours of constant |/?-/? 0 | , a 
determination is made as to how |/?-/? 0 | maps to surface 202. According to the 
formalization, a detector configured to detect scattered electromagnetic radiation is located at 
a fixed point and reflection points for the scattered electromagnetic radiation is moved across 
the surface in such a way as to maintain constant |/?-/? 0 | • For a fixed reflection point, 

25 reflection angles corresponding to constant |/?-/? 0 | trace out a set of observers, such that the 
projected positions are a circle on the surface. That is, the locus of points of constant 
distance to the center is a circle. Since the rules of vector addition do not depend on where 
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vectors are, a locus of reflection points for |/?-/? 0 | will also be a circle (as viewed by a fixed 
observer). 

[0039] The coordinate system, therefore, consists of circles and a radial coordinate, 
wherein the radial coordinate is a function only of |/?-># 0 | . The circles are not necessarily 
5 equally spaced. The center O f of the reflection coordinate system is a point with |/?-/? 0 | = 0. 
The point |/?-/? 0 l = 0 is the specular reflection point. A coordinate system having a center O f 
that is the specular reflection point is referred to herein as the reflection coordinate system 
and is described with the coordinates (r, (f> ). FIG. 6 shows the reflection coordinate system 
(r,<f>) for scattered electromagnetic radiation 215 scattered from an illumination spot 600 that 

10 is illuminated by incoming electromagnetic radiation 205. 

[0040] Another natural choice for a coordinate system is the center O of illumination spot 
600. For convenience illumination spot 600 is shown to be circular, however, this is not 
necessary. The illumination spot may be square or other convenient shape. See, for example, 
the discussion found in Light-Scattering Characteristics of Optical Surfaces, of James Elmer 

15 Harvey, dissertation submitted to the Faculty of the Committee on Optical Sciences, The 
University of Arizona, 1976, herein Harvey '76. The coordinate system centered at the 
illumination spot is referred to that the spot coordinate system. Coordinates for the spot 
coordinate system are written in term of (/?, y/) . The origin O f of the reflection coordinate 

system written in terms of spot coordinates may be written as ( p spec , lf/ spec ) . 

20 [0041] For convenience, the position of a detector 610 (see FIG. 6) is written in Cartesian 

coordinates as ( x 0 , y 0 , z 0 ) . Writing the location of the detector in terms of Cartesian 

coordinates facilitates and simplifies transformations between the reflection and spot 
coordinate systems. The x-axis is set parallel to the line connecting the specular point with 
the detector. The position of the detector is the only point out of the z = 0 plane (plane of the 
25 sample) that is taken into consideration in developing expressions for r and for the reflection 
and spot coordinate systems. For convenience, the angle of incidence of incoming 
electromagnetic radiation 205 and a surface normal, and the angle of reflection of scattered 
electromagnetic radiation 215 and the surface normal is described using the term 0 spec . As 

the angle from the surface normal to the source is known, 0 spcc is taken to be a known 

30 quantity. 
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[0042] Based on the geometry described above, the following relations may be written: 



spec 9 

x Q =p 0 sin y/ 0 , 

*spec = Apec COS Vspec • ^ d 
yO=ys P cc=Po Sin V = PspecS™Vspec- 

From the above expressions, the following equation may be derived: 

/^spec cos ^ = " tan^spec + P 0 cosy/ 0 , and 
A P ecSin^ spec =/? 0 sin^ 0 . 

The following spot coordinates for the specular point may then be derived from the above 
equations and written as: 

pr, sin w n 
tan y/ = — ° 

p 0 cos i// 0 -ztan 9 S 



spec 

Ps PGC =Po -2p 0 z 0 cos if/ Q tan 0^ +z 0 2 tan 2 <9 spec . 



[0043] The limits of the area integration for the BRDF may then be expressed in the spot 
coordinates: 0< r <R where R is the spot radius, and 0< (j> < In . However, the BRDF 
integral separates into reflection coordinates. To simplify the BRDF integral, the limits of the 
area of integration may be expressed in the reflection coordinate system. To write such an 
expression, a transformation between spot coordinates and reflection coordinates is 
developed and is presently described. 

[0044] A given reflection point in reflection coordinates may be written as (r, <f) . The 

Cartesian coordinates for a Cartesian system centered on the specular point are (jc', j/) . In 

spot coordinates, the specular point is (p spec , 0 spec ) = (x spec , >> spec ) , wherein the unprimed 

(x,y) have their origin at the spot center. Then, an arbitrary point in the unprimed Cartesian 
system is 

M a (Wf.) + M 
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In terms of polar coordinates, this may be written as: 

(pcos^,psin y/)=( p spec cos ^ spec + rcos^,p spec sin ^ spec + r sin^) 

Recall that the limits of <f> at a given p , such that p <R , are sought. That is, the limits 

remain inside illumination spot 600. Summing the squares of the x and y components 
5 provides the expression: 

P 1 = +r 2 +p spec r cos ^ cos ^+p spec sin ^ spec r sin^ . 
The condition p < R becomes 



R 2 -p 2 -r 2 

spec - > cos y/ spec cos 0 + sin y/ spec sin 0 . 



[0045] This is a transcendental equation, however, the equation may be solved numerically. 
10 The equation yields the coordinates of the points of intersection between a circle centered on 
the specular point and a circle centered on the illumination spot. There can be zero, one, or 

two such intersections. A parameter /?' = |/? — /? 0 | is defined and an angular parameter a' 
is defined that follows constant /?' contours, a' is chosen to be equal to (f> . Therefore, 
^(r) is the angular length of the arc between the solutions of 

R 2 -p 2 -r 2 

15 = > cos y/ cos 0+sin yr sin^ . 

or 

A^spec 

[0046] While the DSP and BRDF are written in terms of |/?-/? 0 | , the coordinate 
transformations are written in terms of r. By design, the gradient of |/?-/? 0 | is aligned with r, 
so that r may be written in term of r = r( \J3 -J3 0 | ). Furthermore, the geometry suggests that 
the transformation should be monotonic as \J3-J3 0 \ = 0 is the center of the coordinate system. 

20 [0047] The general form of the BRDF 



brdf=— — f f f d ' p ' Hayy* 
P . q hi Ja™ hd dn.dA da d 

in view of the above coordinate transformations may be rewritten as: 
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BRDF= \\*L*P- d p'da', 

JJ A A HO 



6A dQ 

recall that ft' = \j3 — j3 0 \ 9 and a' 15 the angular parameter that follows the constant /?' 
contours. As the integrand is constant along direction J3' , the integrand is parameterized by 
selecting a'=(j). The /? integral may then be written in terms of r, as long as a Jacobian is 
employed: 

BRDF= JJ — — ^@-drd<f> . 
JJ d4 dQ dr Y 

dp' 

[0048] — — tend to be a relatively complicated expression and accordingly a simplification 
dr 

is sought. Without loss of generality, a relationship between r and |/?-/? 0 | is developed for 
<f> =0 as r depends only on |/? — /? 0 | and not on <f> , 



10 \p-P Q \=P'=jsm 2 e+ sin 2 <9 0 -2 sin 6 sin <9 0 . 

[0049] 0 0 may be taken as a known incidence angle and 0 may be written in terms of 
trigonometric functions of the detector position altitude z Q and the (variable) radius r: 



4 + r 2 4z^? 



dp' 

[0050] The derivative may be written as: 

dr 



dP' 



i 



dr 



4 2 z„sinfl 0 . 2 



2, r 2 n r +sin0 o 

K o+r y/r 2 +zZ j 



2 f 2 n ^ 

z 0 r z 0 r sm0 o 

2\(3/2) 



dp 1 dBDRF 

The simplified expression — = for the DSP derived above may now be 

d£l 1J(t)t dr 



wntten as: 
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dp 1 [ dp'\ dBRDF 



dQ IJ(r)r{dr J d/3' 



[0051] FIG. 7 is a high-level flow chart having steps for deconvolving the generalized 
BRDF expression 



according to another embodiment of the present invention. It should be realized that the steps 
shown in FIG. 7 are not limiting on the invention as recited in the claims. Other techniques 
having fewer, substitute, and/or additional steps are within the purview of the invention and 
will be readily apparent to those of skill in the art. Similar to the embodiment outlined in the 
high-level flow chart shown in FIG. 3, at an initial step 700, incident electromagnetic 
radiation is approximated as being essentially collimated. From a physical optics point of 
view, this means that the rays of incident electromagnetic radiation 205 are essentially 
parallel. In view of this approximation, rays of incident electromagnetic radiation emanating 
from a unique differential solid angle dCli are incident on a corresponding unique differential 
surface area dA. FIG. 4 shows rays 405 emanating from a unique differential solid 
angle 410, traveling essentially parallel, and being incident on a differential elements 420 of 
surface 202. As there is an essentially one-to-one relationship between unique differential 
source solid angles dOi and unique differential surface area dA, the contribution to the BRDF 
from these differentials is a delta function for each unique differential source solid angles dCli 
and corresponding unique differential surface area dA. Therefore, the integral over either the 
source solid angle or the integral over the surface area of the sample can be replaced with 
one. In the embodiment presently described, the integral over the source solid angle will be 
set to one. Setting the integral over the differential solid angle d£li to one may be viewed as 
combining the integral of the differential solid angle dQi and the differential surface area dA 
due to the interdependency of these differentials from the collimation assumption. 
[0052] As the integral over the differential solid angle dCli has been set to one, the 
expression for the BRDF has been reduced to a four dimensional integral that may be written 
as 



P. Q . jArea dA d£l d 
At 710, an imaging detector is assumed to collect scattered electromagnetic radiation 215 
scattered from surface 410. An imaging detector is configured to collect radiation scattered 



BRDF=- — — f f f 

P Q JOi JArea kid 
i i 



d% <fe d (Q,,Q rf ,A) 
dQdA dCl d 



dQdAdnj 
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from each surface element dA on the sample surface 202. As light from each surface element 
dA is collected by the imaging detector, collected light may be localized to a very small area 
on the sample. Accordingly, the integral dA over the surface may be estimated as a delta 
function, step 720. That is, the integral over the surface may be replaced with one. The 
remaining integral is over the scattering solid angle dCl d and may be written as: 



BRDF=— — — f 
P. Q,. J ™ 



sin0 d dO d d(/) d 



[0053] A primary aperture of the imaging detector can be described as a set of points 
(0,<f>) in spherical coordinates that are centered on the sample origin O. As (/?— /? 0 ) is a 

function of the set of pointes (9, (f>) , the angles 6 and <p , as well as the constant incident 
10 light angles, a coordinate system may be construct across the primary aperture, such that one 
coordinate follows lines of constant — /? 0 | , while the other is perpendicular to the lines of 

constant \/3 — /? 0 | . A coordinate k x is defined to be the coordinate that follows lines of 

constants \j3 — /? 0 | , and a coordinate k 2 is defined to be the coordinate perpendicular to the 

lines of constant |/? — J3 0 \ . Accordingly, k 2 may simply be written as Ar 2 = |/? — /^o|- 

15 Constant |/? — /? 0 | implies that equation 



/? 0 |=/?' = ^/sin 2 0+ sin 2 0 0 -2 sin#sin 0 O may be written as 
sin 2 6 = 2 sin 0 sin 0 O cos A <f> + C . 
[0054] k x is, therefore, set equal to sin 2 0 . The Jacobian J(k l ,k 2 > ) is calculated to cast 
the BRDF into the coordinates k x and k 2 . The BRDF may then be written in terms k x and 



20 k 2 as: 



brdf=J.*Ma 

*° dQ. v 1 



d(0,<j>) 



d(k x k 2 ) 



dk x dk 2 . 



[0055] D represents the area of the primary aperture of the imaging detector and is 
suitably described in terms of the coordinates k x and k 2 . The Jacobian j(k x ,k 2 ) maybe 
calculated in a customary manner, such as that described in Vector Calculus 4th edition, by J. 
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E. Marsden and A. H. Trombda, published by W H Freeman & Co. April 1996, at pages 372 
376. 

[0056] From the definition of k x , 0 may be written as 0= COS -1 yfk^Afi . <j) maybe 
extracted, resulting in the relationship: 



A^=cos 1 



-k 2 +k x +sin 0„ 



2yfk x sin0 o 



J 



80 _ dd> 
[0057] The derivative — — =0, therefore, — J - need not be calculated. The derivative 



dk, 



dk x 



d0_ 
dk x 



IS 



80 



dk x lyjkQ-k) ' 



and the derivative is 

dk 2 



d(f> 
dk, 



= yj4k x sin 4 0 0 -{k x - k 2 + sin 2 0 0 ) : 



80 8d) 

[0058] The Jacobian is just the product of the derivatives and — — , since the other 



dk x dk : 



80 



term was eliminated by = 0, and may be written as follows: 

8k, 



j{k x k 2 ) = 



d{k x ,k 2 ) 



=(ljk x (l-k x ) (4k x sin 4 0 0 -(k x - k 2 + sin 2 0 O ) 2 ) ] 



[0059] Thus, the expression for the BRDF may be written as: 
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10 



[0060] The limits / and g are those values, for a given k 2y where the contour of |/?-/? 0 | 
intersects the edge of the primary aperture. That is, if (0, 9 $ x ) is the center of the aperture, 
then (#2^2) * s sought, such that 

cos ~ l (sin 0 X cos <f> { sin 0 2 cos <f> 2 + sin 0 X sin <f> x sin 6 2 sin <f> 2 + cos 0 t cos 0 2 ) < a , 

wherein, a is the half-angle of the aperture as viewed from the sample. Expressions for 0 
and <j) in terms of kj and k 2 are given above, and the resulting transcendental equation can be 
solved numerically for a given value of k 2 . Next, the integral over kj can be solved 
numerically, for a given value of k 2 . The BRDF curve is differenced, giving the 
approximation: 



dBRDFjkJ dpjkj 

- ^ K \ K 2), 



dk 2 dd 



wherein /C is the integral over k } 

•k x =g(k 2 ) 



J™ . 
yjk l J(k l9 k 2 )dk l 
k X =f{k 2 ) 



[0061] The following description provides a discussion of a method for implementing the 
above equations for point detectors and imaging detectors to generate a DSP from measured 

15 BRDFs. A method for implementing the above equation for a point detector is described in 
combination with FIG. 8, which shows a high-level flow chart having steps for implementing 
the method according to an embodiment of the present invention. It should be realized that 
the steps shown in FIG. 8 are not limiting on the invention as recited in the claims. Other 
techniques having fewer, substitute, and/or additional steps are within the purview of the 

20 invention and will be readily apparent to those of skill in the art. To implement the described 
method for a point detector, such as a photomultiplier tube, measured BRDF data is 
converted to an equivalent DSP. Input data for the method is measured BRDF data at a set of 
angles, and the output is a calculated DSP at a set of desired angles. 

[0062] At 800, the set of desired DSP angles are selected. At 810, BRDF measurements 
25 are made. The measurements may be made at a specified angular resolution, either by 
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measurement or interpolation. At 820, for each desired DSP angle: (a) a location of the 
specular point is determined, in spot coordinates, making use of equations: 



tan^ s = sin V* 

spec /? 0 cos y/ o -ztsn0 s 



spec 

/4c=Po - 2 Po^os ^ 0 tan ^+z 0 2 tan 2 ^ spec 



which are described in detail above; (b) a radius in specular coordinates for each value of 

|/? — f$ 0 \ that has been measured is determined; and (c) at each of the radii 

[0063] i. dBRDF/djS is approximated by differencing appropriate measured values of 
BRDF, 

dp . 

[0064] ii. — — is calculated making use of equation: 
or 



dr 



zl 2 z sin ft, , 
_Q 5l_ o +sin 2 0 n 



[0065] iii. t(r) is calculated by numerically finding the zeros of equation: 



Psmc T - > cos i// spec cos (f> + sin y/ spec sin <fi , and 



[0066] iv. the DSP at r is calculated making use of equation 



dp 1 



f j o>\ 



dp' 
K dr 



dBRDF 



Each of the above equation in steps (a), (b), and (c) listed above is described in detail in 
above. 

[0067] A method for implementing the above equation for an image detector, such as a 
single pixel of CCD array or CMOS array, is described in combination with FIG. 9, which 
shows a high-level flow chart having steps for implementing the method according to an 
embodiment of the present invention. It should be realized that the steps shown in FIG. 9 are 
not limiting on the invention as recited in the claims. Other techniques having fewer, 
substitute, and/or additional steps are within the purview of the invention and will be readily 
apparent to those of skill in the art. Similar to the method described above for a point 
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detector, measured BRDF data is converted to an equivalent DSP. The input data is the 
measured BRDF at a particular pixel of an image of the sample, at a set of primary aperture 
angles, and the output is a calculated DSP at an arbitrary set of desired angles. 
[0068] At 900, the set of desired DSP angles are selected. At 910, BRDF measurements 
are made. The measurements may be made at a specified angular resolution, either by 
measurement or interpolation. At 920, for each desired DSP angle: (a) calculate the value of 
|/?~/?o| in the cent er of a primary aperture of the image detector; (b) calculate the range of 
|/?-/? 0 | present on the primary aperture, as described by the equations 

tany = P ° Sin ^ 

spec p o cosiy o -ztan0 spec ' ;and 

Apec =Po -2p 0 z 0 cos y/ 0 tan 0^ +z 2 Q tan 2 0 spec 

(c) for each value of |/?-/? 0 | (sufficiently spaced to measure the desired DSP angles): 
[0069] i. approximate dBRDF/d/3 by differencing relevant measured values of BRDF, 
[0070] ii. calculate the limits of the BRDF integral by numerically solving the equation: 

cos _1 (sin 9 X cos <f> x sin 0 2 cos <f> 2 + sin 0, sin cj> x sin 6 2 sin </> 2 + cos 6 X cos 0 2 )<ct 
[0071] iii. calculate K, making use of equation: 



J,k l =g(k 2 ) 
*i=/(* 2 ) 



[0072] iv. solve equation: 



dBRDF(k 2 ) = dp(k 2 ) ^ 



^ \ 2 ) to generate the DSP. 

Each of the above equation in steps (a), (b), and (c) listed above is described in detail in 
above detailed descriptions above. 

[0073] FIG. 10 shows an optical collection system 1000 in which embodiments of the 
invention may be implemented. Optical collection system 1000 may be configured to collect 
BRDF data for calculating a DSP of a sample surface 1010 of a sample 1015. Optical 
collection system 1000 may include an illumination source 1020 configured to illuminate a 
sample with a substantially collimated beam. Illumination source 1020 may be a laser, a 
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mercury vapor lamp, xenon arc lamp, or the like, and may include optical elements 
configured to spread electromagnetic radiation from the illumination source sufficiently 
widely to illuminate all or a portion of sample surface 1010. A detector 1030 may be 
disposed on an extended beam 1040 that is configured to sweep through a variety of angles to 
5 collect scattered electromagnetic radiation at the variety of angles. According to one 

embodiment, detector 1030 is a CCD array or a CMOS array. According to one embodiment, 
extended beam 1040 is between about ten meters and about 30 meters in length, and 
according to a specific embodiment is about 20 meters in length. The extended beam may be 
configured to rotate in angular steps of about 0.05° or more, and according to a specific 

1 0 embodiment may be configured to rotate in about 0.1° steps. 

[0074] A focusing module 1050 may be used to collect scattered light into a given 
reflection angle 6 d to detector 1030. According to one embodiment, focusing module 1050 
includes first and second mirrors 1050a, and 1050b, respectively. Mirrors 1050a and 1050b 
may be parabolic mirrors and may form a Cassegrain telescope, which has a focal plane that 

15 approximately coincides with a collector plane of detector 1030. As the focal plane coincides 
with the collector plane, the formed Cassegrain telescope forms images of the sample surface 
at the image plane. Mirrors 1050a and 105b may be configured as off axis mirrors to 
optimize light directed to the detector. That is, mirror 1050b, which may be configured to 
form a secondary mirror of a telescope, is disposed so as not to block light collected by the 

20 mirror 1050a, which is configured to form a primary mirror of the telescope. According to 
one embodiment, mirror 1050a subtends an angle of about 8.3 X 10" 4 radians or 0.048°. 
According to another embodiment, mirror 1050a subtends an angle of about 2.5 X 10* 2 
. radians or 1.4°. According to one embodiment, detector 1030 is configured to direct 
collected electromagnetic radiation to a computation module 1060 configured to calculate a 

25 DSP from the collected BRDF data. 

[0075] FIG. 1 1 shows an example of a DSP that may be calculated according to 
embodiments of the present invention. Primary mirror 1050a may be positioned in a number 
of overlapping positions and measurements made at each overlapping position in sequence. 
Collected data may be used to achieve a complete deconvolution. The DSP shown in FIG. 1 1 

30 may be generated, for example by interpolating between BRDF measurements. Mismatches 
in the DSP may due to BRDF interpolations made in measured BRDF data. Alternatively, 
the DSP may be calculated from BRDF data that is not interpolated. FIG. 12 shows an 
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example DSP generated from the same BRDF data used to generate the DSP in FIG. 11, 
however, the BRDF data used to calculate the DSP of FIG. 12 was not interpolated. 
[0076] BRDF measurements from which calculated the DSPs (shown in FIGs. 1 1 and 12) 
may be calculated are described below. A substantially collimated beam may be incident on 
5 a sample surface at near-normal incidence (see Fig. 10), such that — fi 0 \ = 0 occurs at 
nearly the same point as 9 = 0 . In such configuration, |/?-/? 0 | is approximated by 

• A 2 " tar S et s P ot subtends an angle of 0.48° (i.e., 8.3 X 10' 3 radians) from a 
detector positioned 20 feet away, for example, on boom 1040. To collect 5 data points within 
0.48° degrees, BRDF measurements may be made at 0.12° increments. For more accuracy, 
10 greater than 5 BRDF measurements may be made, for example, 10 BRDF measurements may 
be may at 0.053° increments. For an incident angle at approximately 45°, near-specular 

measurement density changes by a factor of 1/cos It is noted that these densities of 

BRDF measurements are sufficient for convergence of the deconvolution, however, to get 
resolution on rapid changes of the DSP finer BRDF measures may be made. 

1 5 [0077] Measurement variations (or uncertainties) may constrain the quality of the DSPs 
that are extrapolated from measured BRDFs. Calculated DSP may be affected by 
measurement uncertainty in the scattering angle 0, noise in the power measured by the 
camera, pixilation uncertainty in the solid angle used for the BRDF calculation, uncertainty in 
the out-of-plane angle <fi , and variations in the source power due to fluctuations in source or 

20 differences in the neutral density filters used to filter source radiation. 

[0078] An expression for the DSP uncertainty, which is introduced by uncertainty in the 
scattering angle, is generated via the chain rule by 

9 d(3' 89 

wherein /?'=|/?-/3 0 | (/?' is introduced for convenience of notation). This uncertainty may 

25 be introduced either by uncertainties in the beam positions or by uncertainties in the sample 
angles. 

[0079] For a sample that is a Harvey-Shack material (see Harvey '76) with a power law 
DSP, the DSP is: 

DSP=m(3'-' Y , 

30 such that 
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(sin0-sin0 o cosA<£). 



dDSP _,_ 7 _, 
[0080] By definition, 

\P ~ Po | = ( sin 2 0 + sin 2 5 0 - 2 sin 9 sin 0 0 cos A<f>) ' /2 , 

therefore, 

a/3' = -COS0 

a# " 

[0081] Therefore, the relative imcertainty in the DSP may be written as: 

6 g DSP _ 7m/?' ~ 7-1 cos e (sin 6 — sin# 0 cosA^> ) 

DSP w/?'" 7 " 1 66 

_ 7Cos0(sin0-sin0 o cosA0) 

. ji » 

[0082] In the special case when A0 is zero, |/?-/? 0 | = sin#-sin0 o , the relative 
uncertainty in the DSP may be written as: 

60DSP _ 7 cosfl 5Q 
DSP " sin<9-sin0 o 

[0083] Table 1 shows relative uncertainty in DSP for a power-law index of 7 = 6 and an 
incident angle of 6 0 is about 45°. The angular uncertainty is a function of the incident angle 
at constant \fi - fi 0 1 . The relative uncertainty is within an acceptable margin for 8 = 0.001 . 
According to one embodiment, 8 = 0.001 corresponds to approximately seven pixel widths 
(e.g.; CCD or CMOS pixel widths) of detector 1030. According to one embodiment, pixel 
widths are approximately 13 microns. Table shows the relative uncertainties in the DSP as a 
function of the incident angle and the angular error for such a configuration. 
[0084] In a configuration in which <j> is the vertical displacement 8^ DSP may be 
estimated. That is, 

d/3' _ y/2 sin 0 sin 6 0 sin <f> 
90 ~ J' ' 

such that 
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Table 1 
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[0085] Uncertainty in a CCD array response may affect the DSP measurement in two ways. 
First, inaccurate CCD response can directly alter measured power, and thus change the 
calculated DSP. Second, a first-order uncertainty may be introduced by uncertainties in the 
matching algorithm. If variations in CCD response lead to uncertainties in matching the 
positions of adjacent frames, an effective scattering-angle uncertainty may result. 
[0086] It should also be understood that the examples and embodiments described herein 
are for illustrative purposes only and that various modifications or changes in view thereof 
will be suggested to persons skilled in the art and are to be included within the spirit and 
purview of this application and scope of the appended claims. Therefore, the above 
description should not be taken as limiting the scope of the invention as defined by the 
claims. All publications, patents, and patent applications cited herein are hereby incorporated 
by reference for all purposes in their entirety. 
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